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Abstract 

We reconsider Muskhelishvili-Omnes (MO) dispersive representations of pho- 
ton-photon scattering to two pions, motivated by the very high statistics re- 
sults recently released by the Belle collaboration for charged as well as neutral 
pion pairs and also by recent progress in the determination of the low-energy 
7T7T scattering amplitude. Applicability of this formalism is extended beyond 
1 GeV by taking into account inelasticity due to KK . A modified MO repre- 
sentation is derived which has the advantage that all polynomial ambiguities 
are collected into the subtraction constants and have simple relations to pion 
polarizabilities. It is obtained by treating differently the exactly known QED 
Born term and the other components of the left-hand cut. These components 
are approximated by a sum over resonances. All resonances up to spin two 
and masses up to ~ 1.3 GeV are included. The tensor contributions to the 
left-hand cut are found to be numerically important. We perform fits to the 
data imposing chiral constraints, in particular, using a model independent 
sum rule result on the p e chiral coupling C34. Such theoretical constraints are 
necessary because the experimental errors are dominantly systematic. Re- 
sults on further p G couplings and pion dipole and quadrupole polarizabilities 
are then derived from the fit. The relevance of the new data for distinguish- 
ing between two possible scenarios of isospin breaking in the /o(980) region is 
discussed. 
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1 Introduction: 



Photon-photon scattering into two pions is a process which probes several aspects 
of QCD strong dynamics. In particular, as all participating particles are either 
massless or light, it can probe the low-energy chiral effective theory of QCD. This 
effective theory has now been worked out up to order p 6 [TJ [2] . There are indications 
that at this order, it can represent the exact dynamics to a very high precision in the 
two-flavour expansion. This was shown, for instance, for the tttt 5-wave scattering 
lengths (see e.g. the review j3]). Unfortunately, most of the coupling constants 
of the p 6 chiral Lagrangian are still undetermined. The 77 — > tttt amplitude is of 
particular interest in this respect because of its strong sensitivity to several of these 
p 6 couplings. Physically, these couplings are associated with electric and magnetic 
dipole and quadrupole polarizabilities of the pion. These are important observables 
associated with the structure of the pion. They can be measured, in principle, in 
Primakov experiments or in photoproduction experiments (apart from low-energy 
photon-photon scattering). Such experiments have been performed but the present 
situation is somewhat confused, e.g. the result of MAMI [3] and the preliminary 
result from COMPASS |5J are not in good agreement, to mention only the most 
recent experiments. 

This paper is motivated by the new experimental measurements by the Belle 
collaboration of 77 — > tttt differential cross-sections for both charged pions [HI E] 
and neutral pions [SJ IS]. Only their charged pion results have been used in pre- 
vious theoretical analysis. There has also been very significant progress, recently, 
in measuring the tttt scattering amplitude at low energies by the NA48/2 [101 HI] 
the DIRAC [12] and E865 [13] experiments. We will focus here on relating the 
77 — > tttt experimental results and the low-energy sector of QCD. We will argue 
that using chiral constraints is useful in analyzing the data and that, in return, 
chiral information can be extracted from the data. This might appear puzzling at 
first sight, because Belle's data does not cover the very low energy region: the 7r°7r° 
data covers the range E > 0.6 GeV and the tt + tt~ data the range E > 0.8 GeV. 
Extrapolation is possible due to theoretical properties of scattering amplitudes in 
the standard model, in particular, the property of analyticity of partial-wave am- 
plitudes as a function of energy. This property, we recall, is a proved consequence 
of confinement in QCD [13] . Combining with unitarity of the S'-matrix enables one 
to disentangle the effects of the final-state interaction by the Muskhelishvili-Omnes 
(MO) method [151 EE]- Application to 77 — > tttt amplitudes was discussed for the 
first time in ref. [T7]. Explicit results for 77 — > tt + tt~ taking into account current 
algebra constraints were obtained in [T8"| H9] . This was reconsidered a few years 
later [201 Ell 122] after the first reliable experimental results in the low-energy region 
became available. Refs. [211 [22] also discuss how the MO dispersive representa- 
tion matches with the chiral one-loop representation, which had been computed in 
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refs. [231 [21] and eventually lead to a parameter free prediction in the low-energy 
region. 

A simplifying feature of low energy is that tttc scattering can be considered elas- 
tic. If one is interested in the 1 GeV region or slightly above, it becomes necessary to 
take inelasticity into account. This is feasible due to a specific feature of tttt scatter- 
ing: the fact that inelastic scattering to Att or 6tt states (which are not treatable by 
MO methods) are suppressed in practice and can be neglected up to E ~ 1.2 — 1.3 
GeV. The remaining relevant inelastic channels, KK or rji], are two-body channels 
which are perfectly treatable, in principle, by MO methods. In this paper, we take 
into account — > KK scattering, which is particularly important in the 1 = 
S-wave near 1 GeV. This amplitude suffers from a long lasting unresolved experi- 
mental discrepancy very near the KK threshold which limits the determination of 
the properties of the / (980) scalar meson. It is interesting that this discrepancy 
can be interpreted in terms of two different scenarios for isospin breaking. We will 
discuss the relevance of Belle's results near 1 GeV in eventually clarifying this issue. 

Application of the multichannel MO method to 77 — > tttt has been attempted 
first in ref. |25J and then discussed in some detail, on the basis of specific T-matrix 
parametrizations for which the MO integral equations are analytically solvable, in 
ref. [26J. More recently, it was applied to Belle's data on 7r + 7r~ in ref. [27J. These 
authors have considered the extrapolation of the amplitude in the complex plane, 
so as to define, and then extract, the couplings of the scalar mesons a (600) and 
/o(980) to two photons. Determination of these couplings had aroused significant 
interest in the literature (e.g. [231 1221 EDI EI], a more complete list can be found 
in [27]). One motivation is to probe the structure of scalar mesons and identify the 
glueballs. There is some scatter in the results obtained. In the present paper, we 
concentrate on extrapolating on the real axis, towards the low-energy region. This 
brings constraints on the amplitude which should prove useful also for extrapolating 
away from the real axis. This will be discussed elsewhere. 

The plan of the paper is as follows. After introducing notation for the amplitudes 
and their partial-wave expansions we write the unitarity equations in the one and 
two-channel approximations. Next, we formulate the MO-type dispersive represen- 
tations. Concerning the left-hand cut, we find it advantageous to treat differently 
the QED Born term and the multipion contributions in the MO representations. 
The latter are kept in the form of a subtracted left-cut spectral integral. We then 
implement the (usual) approximation of retaining only resonance contributions. All 
resonances with mass up to ~ 1.3 GeV are included and we show that a certain reg- 
ulation operates between resonances of different spin and different parity depending 
on the helicity states. Subtractions at s = are introduced in the MO dispersive 
representations in order to suppress higher energy regions in the integrands where 
our truncated unitarity equations no longer apply. The subtraction constants have 
simple relations to dipole and quadrupole pion polarizabilities and are to be deter- 
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Figure 1 : Mandelstam plane showing the physical regions for 77 — > 2ir and '-fir — > 77r 
processes. 

mined from fits to the data. Chiral constraints may be applied to the fit. This 
is necessary because the errors in Belle's data are completely dominated by sys- 
tematics and the usual statistical interpretation of the \ 2 does not apply, strictly 
speaking. Only one of the relevant p 6 chiral coupling constants is known in a model 
independent way from a chiral sum rule. We show that this information implies a 
relation between dipole and quadrupole polarizabilities of the neutral pion which we 
implement in the fit. Then, we describe our inputs for the tttt — > tttt, KK T- matrix 
elements. Our subtracted dispersive integrals emphasize the low-energy part of the 
integrands. We employ a parametrization which allow for some freedom near the 
KK threshold since the Belle data probes this region in some detail. Finally, we 
display comparisons between the fitted MO amplitudes and the experimental data 
and discuss the implication for the pion polarizabilities and the p 6 chiral coupling 
constants. 

2 Kinematics, unitarity relations 

We consider the processes 7(<?i, A)7(g 2 , A') — > vr + (p 1 )7r _ (p 2 ) or vr ^!)^ ^) where 
A(A') = ±1 are the photons helicities. We will also consider j(qi, A)7(g 2 , A') — >• 
K(pi)K(p2) with 1 = which plays a role in the S"-wave via coupled channel uni- 
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tarity. We take the Mandelstam invariants as 



s = (qi + q 2 ) , t = (q 1 -p 1 ), u = (gi - p 2 ) (1) 

which satisfy s + t + u = 2m 2 . The physical regions in the Mandelstam plane for 
77 — > 2tt and for the crossed-channel amplitude 77? — > ^ix are shown in fig. [TJ The 
scattering angle 9 in the two photon center-of mass system is related as follows to 

s, t, u 

t — u 

cos 6 = : . (2) 

's(s — Ami) 



We write the S"-matrix element for 77 — > it + tt as 

< , ^( Pl )n-(p 2 M qi , A) 7 (g 2 , A') > = ie 2 {2v)H\P f - H c xy (s, t) (3) 



out 



factoring out e 2 as well as the explicit dependence on the azimuthal angle (p. With 
this convention, H\\i is a function of the Mandelstam variables. Similarly, in the 
case of 77 — > 7T 7r° we denote the scattering amplitude by H xx ,. In the case of 77 — > 
KK scattering we denote the charged and neutral amplitudes by K xx , and K xx , 
respectively. We assume that isospin is exactly conserved by the strong interaction. 
It is then useful to consider the amplitudes which correspond to ttit or KK final 
states with definite isospin I. We will label them as H xx , and K xx , . Because of 
parity conservation, only final states with even values of the angular momentum J 
are allowed, J = 0,2,4, ■ • • Invoking also charge conjugation invariance the isospin 
values must be / = or / = 2 in the case of 7T7t, while in the case of KK both 1 = 
and 1 = 1 can couple to 77. The relations between the amplitudes 77 — > 7r + 7r~, 
7r°7r° and the isospin ones 77 — > (7r7r)/ =0) 2 read 
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In the case of kaons, the analogous relations read, 




(4) 



K c X x> ) I v 2 
K% I /X 




(5) 



It is useful to carry out a tensorial decomposition of the photon-photon scattering 
amplitudes. Writing 

H xy ( qi , Pi ) = e?(AK(A')HVfe, gi ) (6) 
where are the polarization vectors of the photons. W^ v can be decomposed as 

= A(s, t, u) 2V + B(s, t, u) T 2ilv (7) 
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where 



lflU 



-, s 9iiv — llv&fi 



T 2ilv = 2s -{t- ufg^ -2{t-u) {q lv ^ - q 2 ^ v ) (8) 

where A = (pi — p 2 ). In this manner, the Ward identities are satisfied as follows, 

(tfWV = <&Wp, = . (9) 

The functions A and B satisfy analyticity properties as a function of s, t, u and 
they are symmetric under crossing (t, u) — > (u, t) (because of Bose symmetry of the 
two photon system). One can express the helicity amplitudes in terms of A and B 
as 

followfl 

1 

( — l-l 



II = l-l A- s(s - Ami) B 



H + _ = H_ + = 4{tu -m*)B . (10) 
Finally, the differential cross section for 77 — > mi has the following expression 



da 



a 
8s 



f3 n (s) (\H ++ \ 2 + \H + _\ 2 ) 



with 



1- 



Ami 



(11) 
(12) 



2.1 Partial- wave expansions 

In order to perform the partial-wave expansion for helicity amplitudes we use the 
Jacob and Wick [33] formulas, 

(8(pX c \d\T\0,0, \ a \ b ) = iV^(2J + l)D*/ a _ Xb Xc _ Xd (9, (p)(\ c \ d \Tj\\ a \ b ) (13) 

j 

where N is a normalization factor which can be chosen arbitrarily. Let us list below 
the partial-wave expansions for all the scattering amplitudes which are relevant in 
our work, 



77 — 7- 7T7T : 


H\\< 


77 — > KK : 


K{x> 


7T7T — > 7T7T : 


F 1 


7T7T — > KK : 


G 1 


KK KK : 


R° 



E(2J+1) h^^dty^O) 

^E(2J+1) k^y^dty^O) 

327TE(2J+l)/i( S )rf J W 



16ttE(2J + 1) r°j(s) d J oo (0) ■ 



(14) 



lr The polarization vectors are chosen in accordance with the phase convention of Edmonds 
for spherical tensors. 
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The different normalization factors are chosen such as to ensure simple formulas 
for the unitarity relations satisfied by the partial-wave amplitudes. The general 
unitarity relation reads, for a given T-matrix element 

T fi ~ T if = i Yl T nf T ni ■ (15) 
n 

If the energy is sufficiently small, the sum over intermediate states is limited to 
just one state, im (elastic unitarity). We will consider this to be a reasonably good 
approximation for our purposes except in the I = J = case (see sec. At the 
level of the partial-waves, if J ^ or I ^ we then have 

Imh^is) = 0(s - 4m^( s )/j( s )/^ A ,( s ) ■ (16) 

For I = J = we also include KK in the sum (j!5p . The unitarity relation can be 
written in matrix form 

/ Imfc° ++ ( S ) \ _ / h$ ++ (s) \ 

\ Imfc° ++ (s) J V K°++( s ) ) { V 



with 



T _(f8(s) 9°o(s)\ y _ (0(s - AmDMs) \ 

" I 9°o(s) r °( S ) h -[ 9(8 - 4ml)P K ( S ) i ' ^ 



3 (Modified) Omnes -Muskhelishvili representa- 
tions 

The partial- wave photon-photon scattering amplitudes hj XX ,(s) are analytic func- 
tions of the variable s with two cuts on the real axis: 1) a right-hand cut extending 
from 4.m\ to oo and 2) a left-hand cut extending from — oo to 0. The discontinu- 
ity along the right-hand cut is given by the unitarity relations (TIB]). ffT7|) . These 
properties are the basis of the Muskhelishvili-Omnes (MO) method for treating the 
final-state interaction problem, which has been applied for the first time to the 
77 7T7T scattering amplitudes by Gourdin and Martin [TTJ . The usual method [TS] 
is based on writing a dispersion relation for the function, 

F(s) = n-\s)[F(s)-F L (s)} (19) 

where F is the amplitude of interest, Fl the part of this amplitude which has a 
left-hand cut and Q is the Omnes function. The function F(s), by construction, has 
only a right-hand cut. We will use a slightly modified version here, which treats 
on a different footing the part of the left-hand cut associated with the QED Born 
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term (which is known exactly) and the remaining part, which we will associate with 
resonance exchanges. Instead of fTiTJl) . we will consider 



(8)[F{8) - F*™{8)] 



(20) 



i.e. we subtract only the Born term piece such that the function F mo d(s) has both a 
right-hand cut and a left-hand cut. A simplification arising in this modified approach 
is that all polynomial terms are absorbed into the subtraction constants of the 
dispersive representation. 

In practice, in order to evaluate the imaginary part of F(s') (or F mo d{s')) on the 
right-hand cut, which is needed in the dispersion relation, one sets Im [fi _1 (s')F(s')] 
equal to zero in eqs. (1191) (1201) . This is exact in the energy region where scattering 
is elastic but becomes inaccurate at higher energies. The influence of this inaccu- 
racy can be reduced by writing down over-subtracted dispersion relations, so as to 
suppress the integrand in the inelastic region. With this in mind, it also makes 
sense to define the Omnes function over an infinite range, making a plausible guess^l 
concerning the behaviour of the phase-shift as s' — > oo. Alternatively, one could 
perform all the dispersive integrations over a finite range s' < s c and approximate 
the contributions from the range s' > s c by polynomials with unknown coefficients. 
This should be practically equivalent to the procedure adopted here. In the end, the 
sensitivity of the results on the higher energy ranges of the various integrals would 
have to be included in the errors. 

Let us now consider the MO representations in more detail. 

• S-wave 1 = 

For the 1 = S-wave, it is necessary to generalize the MO representation to two 
channels in order to properly describe the /o(980) resonance energy region. The 
Omnes function must be replaced by a 2 x 2 Omnes matrix 



The matrix elements of the Omnes matrix are analytic functions of s with only a 
right-hand cut as in the one-channel case. The discontinuities along this cut are 
given, in terms of the 2x2 T-matrix, by equations analogous to (ITT)) that read, 



Unlike the one-channel case, the MO equations have no known analytic solutions for 
two or more channels |16j . but accurate numerical solutions can be constructed [Ml 

2 To be more specific, we assume that 1 = 2 phase-shifts tend to to zero. For 7 = and J = 
we take: lim^oo S n „(s) = 2n, Iim^oo 5 K k{s) = and lim^oo \T nw ^. K ^(s)\ = while for I = 
and J = 2 we assumed that the phase-shift goes to n. 




(21) 




(22) 
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The MO representation couples 77 — >■ (tttt) i=0 S-wave amplitude and the 
77 — > (KK)j =0 S-wave amplitude. We write, a priori, a representation which 
involves four subtraction parameters 

0,Born,_x \ \fb^S + b'^S 2 

+ Q(s) x 



h° 0t++ (s) \ _ / ^ 



bPs + bTs 2 




s 3 f~ s o ds' =-1 

+ W-oc {S'f(s'-S) 

-- f 2 7 .3^ J mH'^O ( j^trn) S f { ) ■ (23) 

As usual, subtraction constants have been set equal to zero in order to comply with 
the soft-photon theorem [36j El] near s = e.g. 

&O.++00- Wto = ( 24 ) 

The left-cut functions foo++ n (s'), ^o'++ n ( s ')j ^o!++( s ')5 ^o,++( s ') which enter into 
this representation will be given explicitly in sec. SJ 

Due to the property of real analyticity the discontinuity of the amplitude ++ (s) 
across the right-hand is expressed in terms of a phase (poo 

h° , ++ (s + *0 = e 2i ^h° 0)++ (s - ie) (s > Ami) • (25) 

This phase is equal (modulo it) to the inr phase-shift below the KK threshold by 
Watson's theorem. The representation given above provides a modelling of 0oo above 
the KK threshold (depending on the polynomial parameters) which is plausible 
below the effective onset of 47r inelasticity. The amplitude ft-o ++ (s) satisfies a one- 
channel Omnes representation in terms of the Omnes function associated with </>oo 
and two polynomial parameters. We have verified this property as a check of our 
numerical calculations. 

• S- wave 1 = 2 

In this case, KK inelasticity is not allowed. We will then disregard inelasticity in 
the energy region of interest. We write an MO representation with two subtraction 
constants 



hl ++ (s) = hl'++ n (s) + Q 2 Q (s) 



b^s + b'^s 2 (26) 



s 

+ ~n J-00 ng(s')(s') 3 (s'-s)""' ' ir Jiml \n 2 (s')\(s') 

Here 5q(s) is the I = 2 tctc phase-shift and the Omnes function is given in terms of 
Sl(s) by ' 

a ^^Cw^ d3 ')- (27) 
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• D-waves 1 = and 1 = 2 

In the 1 = case, the partial-wave analysis performed by Hyams et al. [3H] found 
the inelasticity around the /2(1270) peak to be of the order of 30%. The PDG [39J 
now quotes a smaller value, approximately 15%, of which only 5% is due to KK, the 
remaining part being due to An. We will therefore not attempt a coupled channel 
description in this case and essentially ignore the inelasticity. In the 1 = 2 case, the 
final-state interaction is very small and will also ignore the inelasticity. The Omnes 
method differs from that discussed above only by the fact that we must account 
properly for centrifugal barrier factors. At small energies, indeed, it is not difficult 
to see that the amplitudes with J = 2 should behave as follows, 



i I, Born/ \ 

i I, Born/ \ 
-«2+- 0) 



s 2 (s — Ami] 
Ami 



sis 



{21 



This is implemented by multiplying the Omnes function by (s — Ami) and by setting 
the appropriate subtraction constants to zero in the dispersive representation. One 
then obtains, 



Ami 



x 



,(-0 



+ 



-SO 



Im h 



1, Res i 

2, ++ \ 



nJ-oo ^(s')(s') 3 ( s/ - 4m£)(s' 



-ds' 



s 

+ - 



sin dl(s') h 2 



I, Born/ /\ 
S 



71 J Ami |fi£(s')|(s') 3 (s' - K)( S ' - S ) 



ds' 



(29) 



and 



hi 



hi>* ™(s) + nl(s)s(s-Aml 
~ s o Im h 2 



d^ + 



I,Res (sl) 



S '" 00 
+ - 



7r7-oo ^(S')(S') 2 (S'-K)( S '-S) 



ds' 



sm o [s ) ti 2+ - ( 



ds' 



(30) 



These representations involve four subtraction constants and S 1 '. 

Let us finally remark that polynomials can be introduced in MO representations 
in different ways (for instance multiplying the Omnes functions). It is easy to show 
that any representation can always be recast in the form given above. These are 
convenient because the relations between the polynomial coefficients and the pion 
polarizabilities are now particularly simple. 
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3.1 Sum-rules 



Above, we have written over-subtracted dispersion relations in order to suppress 
the contributions from large values of \s'\ in the integrands. Such representations 
can be valid only in a finite energy domain since, in general, they lead to diverging 
amplitudes^ when s — > oo. A priori, the exact asymptotic behaviour of the partial- 
wave amplitudes is not known, but S'-matrix unitarity provides the following bound 



IM 



J.XX' 



< 



IQtt 



(31) 



which implies that one could write representations with fewer subtraction constants 
if one assumes that the integrands are known sufficiently well. This amounts to 
expressing some of the subtraction constants as sum rules. Such sum rules are very 
simply obtained from the representations written above by requiring that the most 
singular terms as s — > oo are cancelled. For example, the .D-wave constants 
get expressed as 



SR 



-(I) 
J 3,++ 



(I) 



SR 



with 



00 

J n,XX' 



lmh 2.XX'( S ') 



R 



(i) 

n,XX' 



nJ-oo ^(s')(s') n (s'-K 
sin Sj(s') hj^js') 
n hm\ |^(s / )|(s , ) n (s' -Ami 



-ds' 



(32) 



-ds' . 



(33) 



A useful test (which we will perform) of the validity of the integrands, in particular 
of the modelling of the left-hand cut, is to verify that the values of the subtraction 
constants obtained from such sum rules are not significantly different from those 
obtained from fitting to the experimental data. 



4 Left-hand cut 

In order to proceed with the previous formulas we need to specify the left-hand 
cut pieces of the 77 — > tctt amplitudes (as well as the analogous contributions to 
77 —> KK in the case 1 = 0, J = 0). Quite generally, the left-hand cut can be 
associated with singularities in the crossed-channel (i.e. 77? — > 771") partial-waves. A 
derivation based on the Mandelstam double-spectral representation can be found in 

3 The asymptotic conditions implemented in the T matrices imply that the 1 = and 1 = 2 
Omncs functions behave respectively as s -1 and s° at large s. 
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ref. [17J. The first cross-channel singularity is the pion pole, followed by the unitarity 
cuts due to 2tt, 3tt etc... Here, the pion pole contribution will be exactly taken into 
account. The discontinuities p n7r (s) associated with the unitarity cuts are calculable 
from ChPT at small s, but they are strongly suppressed in this region: p^ w has chiral 
order p 6 while p2 n is even more suppressed and has chiral order p s . Here, due to 
the lack of detailed experimental information on the 77r — » 771 partial-waves, we will 
content ourselves with simple resonance approximations. 



4.1 QED Born amplitudes 



Let us first consider the QED Born term contribution to 77 — > tt + tt (and 77 — > 
K + K~). The standard result for the helicity amplitudes reads 



H 



c,Born 
-+ 



2smz 



(t - ml){u - ml) 



H 



c,Born 



2(tu - ml) 



(t-ml)(u-ml) 



(34) 



The corresponding Born amplitudes for 77 — > K + K~ are, of course, the same, with 
m n replaced by mx- Performing the partial- wave projection according to eqs. f |T4l . 
one obtains for J = and J = 2, 



ic,Born/ \ 

K++ 0) 



1 c,Born/ \ 



j c,Born/ \ 

hi, is) 



S 

2ml 



#00 



LJs) 



ms) 



with 



V6 
4 



LJs) 



ms) 



ti(s)L w (s) 



2 10 

W) + ^ 



1 . 1 + Ar(s) 

log 



(35) 



(36) 



We will also need the J = partial-wave amplitude for 77 — > K + K~ which, as a 
result of the normalization (j!4p reads 



ic,Born/\ 



4y/2m 2 K 



L K (s) 



(37) 



and the isospin 1 = and 1 = 2 projections which are easily deduced from (j2J), (jSJ). 



4.2 Resonance contributions 

Let us review below the modelling of the left-hand cut as a sum of resonance pole 
contributions. Resonances which can contribute must have spin larger than or equal 
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to 1. We will consider vector, axial- vector, tensor and axial-tensor contributions. 
We start by determining the form of the amplitudes and the relation of the coupling 
constants to the radiative decay widths. We will then discuss the phenomenological 
determination of the coupling constants. 

• Vector resonances 

We can start with a Lagrangian coupling a vector meson V a pseudo-scalar meson 
P (a pion or a Kaon) and a photon field of the form, 



C VPl = eC v e^F^d a PVp . 
After a small calculation we reproduce the result first obtained by Ko 



c 



v 



v 



m?/ — t 



(s - Am 2 P - 4t)TT + ~7T 



+ u), C\ 



1 

-i 

2 



-r 2 



(3? 



(39) 



The relation between the coupling Cy and the decay width of the resonance reads 

(m v — m 2 P ) 3 



r 



3m? 7 



(40) 



Performing the partial-wave projections of these amplitudes one finds for J = and 
J = 2 



AG 



v 



m\ 



p n (s 

1 



(s) 



C V ^-sPk(s) 



-L v (s) + s 

-3X v (s))L v (s) + 6X v (s) 



;i - X 2 V (s)f L v (s) + -X v (s)(5 - 3X v (s)) 



(41) 



where the logarithmic function Ly reads 

Xy(s) + 1 



Ly(s) = log 



Xy(s) 



Xy(s) 



2m 2 v 



2m 2 + s 



sfiJs) 



(42) 



One remarks here that a term linear in s appears in the expression for the amplitude 
ho ++ . This is an illustration of possible polynomial ambiguities. Indeed, if one uses 
an antisymmetric tensor description of a vector particle (e.g. [41] ) this term would 
appear with the opposite sign. In the modified MO representation we need only the 
imaginary parts of the partial- wave amplitudes along the cut, which are free of any 
ambiguity. The left-hand cut is contained in the function Ly(s). Rewriting Xy as 



X V (S) 



1 + 



Am v (s + s v ) 
s(s — Am 2 ) 
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[m v — m 

9 

mi 



2\2 



(43) 



one deduces that the cut extends from — oo to — sy. One then derives the imaginary 
parts of the amplitudes 



7T 



— ImhY i 

7T A_l 



-Im/ji 7 , ( 

7T A+ 



-4a 



-2Ct 



-s - Sy, 



-^Xl{s)-\)Q{-s-s v ) 



s v ) 



(44) 



We remark that the sharp cutoff 9{—s — sy) which appears in the imaginary parts 
is associated with the narrow width approximation. Introducing a finite width 
smoothes the cutoff in the region s = —sy. The corresponding contributions to 77 — > 
KK partial-wave amplitudes (i.e. ko ++ , k^ ++ ko + _) are obtained by replacing in 
eqs. (jUJ m n by m K and Cy by y/2C v (in accordance with the normalizations defined 
in eqs. f TT4|) ). 

•Axial-vector resonances 

We will designate generically the C-odd axial-vectors as B and the C-even ones by 
A. The following Lagrangian describes B — > P7 coupling, 



/W 7 = eC B F^B v P 



(45) 



from which we easily deduce the amplitude 77 — > PP amplitude corresponding to 
B exchange, 



c 



B 



B 



m\ — t 



[s - Am 2 P + 4t)?T + ~Tf 



+ (t^u) , Ci 



-Cl 



with 



(46) 



(47) 



As before, the couplings C B must be deduced from experimental data on B — > P7 
decay using 

r (m| -mp) 3 

1 B^-P-y 



aC 



B~ 



3m 3 B 



(4f 



Working out the helicity amplitudes corresponding to eq. (|46|) . one sees that they 
are simply related to the helicity amplitudes associated with vector resonances: 



H^ + (C B ,m B , s,t) 



-H+ + (Cy — »■ Cb , — > m B , s, t) 
H+_(C V -» C*b , my -»■ m B ,s,t) 



(49) 



from which one easily deduces the imaginary parts of the partial-wave projections 
from f|44|) . 
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• Tensor resonances 

The Lagrangian coupling a tensor meson (quantum numbers 2 ++ ) to a photon and 
a pseudo- scalar meson must have the following form, 

C TPl = eC T e^F^T a x d x dpP . (50) 

From this, we can first deduce the relation between the decay width and the coupling 
constant 



r ^ = ^ V " J ' > =7?- ( 51 ) 



[m l - mlf ~ = q 
hm\ ' T 16 



After a small calculation we can obtain the form of the 77 — > PP amplitude gener- 
ated by tensor meson exchange 



Wf = A T (s, t)T? v + B T (s, t)Tf + (t «-> u) (52) 



with 



T{,) 1m\{m\-t) 

A T (s, t) = 2(s - Ami - 4t)B T (s, t) - 8 ^ ~<) 2 (53) 

[m A T — t) 

From here, we can construct the helicity amplitudes and their projections. We quote 
below the imaginary parts of these along the left-hand cut 

I — ~ 771 2 

-Im/^ ++ (s) = -4Ct^ [4s + 3s T )] 9(-s - s T ) (54) 
-Im/^ ++ (s) = -2C T ^f-(3X^(s) - 1) [4s + 3s T ] 9(-s - s T ) 
-Im h T 2+ _{s) = y^C T s^(s)(X 2 T (s) - l) 2 [4s + s T )} 9(-s - s T ) . 

71 4 

• Axial-tensor resonances 

For illustrative purposes, finally, let us consider axial-tensor (i.e. with quantum 
number J p = 2~) resonances. The relevant Lagrangian has the following form 

C TaPi = eC TA F^Tfd^P (55) 

and we deduce the following relation with the radiative decay width 



{m\ A -m P f - _C\ A 
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Next, one computes the diagrams contributing to 77 — > PP and one finds the result 



G 



(t + ml) 2 + Am^ A (s - 



m„ 



A TA (s,t) = 2(s-4ml + At)B TA {s,t) 



8C TA (t-ml) 
t) 



(57) 



Constructing the helicity amplitudes, one notices the following simple relations be- 
tween the tensor and the axial-tensor amplitudes 

H1\{C Ta , m TA ,s,t) = -Hl + (C T C Ta , m T ->■ m TA ,s,t) 

Hl± (C Ta ,m TA ,s,t) = Hl_ (C T -> C Ta ,m T ^m TA ,s,t) (58) 

Alternative expressions can be derived for the resonance exchanges, which are some- 
what more transparent physically and easier to generalize to arbitrary spin S. They 
involve rotation functions df y( z R) where zr is the center-of-mass scattering angle 
for the crossed- channel amplitude 77? — > 77? at the resonance mass, 



1 + 



2st 



(t - ml) 



t=m 2 R 



1 + * 

sr. 



(59) 



(using (I43p ) and generalized Legendre polynomials PJ 1 (Xr). The vector-exchange 
partial-wave helicity amplitudes can be expressed as 



lmh v 



0,++ 



» 



-4CV, 



mySy 



— IrahX 1 
71 AH 



[s) = -AC X 



s(s — Am\ 
myS V 



-.P^X v )dl_ 1 (z v )6(-s-s v ) 



S(S 



4ml 



--P%{X v )d\Az v )e{-s-s v ) 



71 



VQ^ rn 2 v s v 

• >s(s — Ami' 



P 2 (X v )d 1 ll (z v )e(-s-s v ) . (60) 



Analogously, the tensor exchange amplitudes involve the rotation functions d\ y{zr) 



2 c ,2 



71 



Im h o,++( s ) 



-Im h 2>++ (s) 



-AC 7 



-4C 7 



m^rS 



s(s — Ami 



=P°(X r K_ 1 (z T )0(-s-s T ) 



S{ S 



Ami 



2 c .2 



-.p^x^dl^iz^ei-s-ST) 



-Im hi 



mrpS 



71 



^ - 1 s(s — Ami 



P%{X T )dUz T )e{-s-s T ) • (61) 
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From these expressions, it is easy to guess the form of the amplitudes for an arbitrary 
angular momentum J, or generated by the exchange of a meson of arbitrary spin S. 
Because of the polynomial functions df x ,(1+2s/sr), the behaviour of the amplitudes 
as a function of s becomes worse as the spin S of the exchanged resonance increases. 
As usual, one expects that a Regge-type regularization will occur upon including an 
infinite set of resonances. In practice, one can simulate this by introducing a cutoff 
on the left-cut integration. 



AA' 


VAT Ta 


++ 
+- 


- + + - 

- - + + 



Table 1: Signs of the imaginary parts of the helicity amplitudes as generated from 
various resonance exchanges 

Let us make a remark on the signs of these resonance amplitudes. Varying s 
from — oo to —sr, the quantity Xr(s) varies from —1 to 1. The Legendre polyno- 
mial P2(Xfj) passes through a zero, implying that Im/i^ ++ changes sign while the 
amplitudes Im/i^ + _ and lmh,Q ++ do not. This feature partly explains why, upon 
integration over s, the H — P-wave amplitude is larger than the ++ one. Further- 
more, there are alternating signs between the various resonance contributions: see 
table [TJ The table shows that for the ++ amplitudes the signs alternate between 
resonances of different spin and between P-even and P-odd resonances while for he- 
licity H — amplitudes the signs alternate between resonances of different spin. The 
resulting behaviour of some of the integrands from which the left-cut functions are 
computed is illustrated in fig. [2J This figure illustrates the numerical importance of 
the tensor contribution in the helicity H — amplitude. 

4.3 Phenomenological determination of the coupling con- 
stants: 

• Neutral resonances 

The neutral resonances which can decay into a photon and a ir° must be odd under 
charge conjugation. This is the case of the vector mesons p°, u and their properties 
are rather well known experimentally. The results from the PDG [39] and the 
corresponding values of the couplings Cy are collected in table [2] below including 
also the result for the if* (892). 

The experimental information concerning the C-odd axial meson radiative de- 
cays is not as detailed as in the case of the vector mesons. The PDG quotes a result 
for 6i"(1235) decay: T(bf (1235) — > 7i +/ ~f) = 240 ± 60 KeV, whereas the correspond- 
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-4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 -10 -8 -6 -4 -2 

s" (GeV 2 ) s' (GeV 2 ) 



Figure 2: Left-cut integrand for ++ (0) (left figure, see eq. ( T23]) ) and h® + _ (0) (rigM 
figure, see eq. / f30|) J illustrating the role of the various contributing resonances. 



ing radiative widths of the neutral axials 65(1235) and /ii(1170) have not yet been 
measured. A rough estimate of these, using nonet symmetry, is: 

C&0(1235) ^ C&+(1235)> C h o (1170) ~ 9C b + (1235) . (62) 

Concerning the strange axials, we can use here the experimental results on the 
radiative decays of the neutral i^i(1270), f^i(1400) which can be found in the PDG. 
We collect this information and the results for the Cb couplings in the lower part of 





T (KeV) C fi (GeV" 2 ) 


bJ 7T°7 
P° -)■ 7T°7 
(f) — > 7T°7 

K*° K°~f 


703 ± 25 0.66 ± 0.023 
89 ± 12 0.09 ±0.01 
5.4 ±0.5 (0.2±0.02)10~ 2 
117 ± 10 0.20 ±0.02 


/ii(1170) ->• vr°7 
foi(1235) ->• vr°7 
Ki(1270) ->• K°7 
Xi(1400) ->• K°7 


~ 0.45 

~ 0.05 
73 ±29 0.024 ±0.010 
280 ±46 0.063 ±0.010 



Table 2: Radiative widths of neutral vector mesons and of neutral C—odd axial- 
vector mesons from the PDG PJ^|/ and the corresponding couplings Cr. 
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table [2J Finally, both C-odd and C-even axial-tensor mesons are expected to exist 
in the quark model with a mass around 1.7 GeV [12]. Experimentally, the C-even 
axial-tensor meson 7^(1670) is mentioned in the PDG, but not the C-odd. The 
decay width of the 7^(1670) into 77r is not known. 

• Charged resonances 

Here, we can have contributions from charged vector mesons p + (770), iC + (892)) and 
from charged C-odd and C-even axial- vector mesons ct+(1260), 6^(1235) ,7^(1270), 
i^i~(1400). We have also considered the contributions from charged tensor mesons 
ct2~(1320), i^2 + (1430) since their masses are comparable to those of the axial- vectors. 
The relevant couplings Cy for the vector mesons in the charged case can be deduced 
from experiment and are collected in table [31 In the case of the axial- vectors, results 





T (KeV) C R (GeV- 2 ) 


p + — > 7T + 7 
K *+ R + 1 


68 ± 7 0.066 ±0.007 
50 ± 5 0.085 ± 0.009 


0^(1260) tt+7 
67/(1235) -> tt+7 
^+(1270) fC+7 
^+(1400) K + ~f 


640 ± 240 0.15 ±0.06 
230 ±60 0.05 ±0.01 
~ 0.20 
~ 0.00 


4(1320) -> tt+7 
K 2 *+(1430) -> K + 1 


287 ±30 0.052 ±0.005 
241 ± 50 0.053 ±0.011 



Table 3: Same as table for charged vectors, axial-vectors and tensor resonances. 

are available for the &i~(1235) [43J as well as for the (1260) jH] from Primakoff 
experiments. We must however keep in mind that results from photoproduction 
experiments [HI [46] suggest that the radiative width of the af (1260) could actually 
be smaller than claimed in [H]. Concerning the charged strange axials, i^i(1270), 
i^i(1400), unfortunately, no experimental information is available on their radiative 
widths. Rough estimates can again be made using nonet symmetry, which leads to 
the following relations 

CV<>(i270) ± CkO(i4oo) = 4C b + (1235) 

^V+(1270) + ^+(1400) = ^6+ (1235) + ^a+(1260) ' (^) 

The first relation is obeyed by the experimental results within a factor of two. In 
order to determine the couplings C K +^ 1270 ^ and C K +^ U00 ^ separately we note that 
one of them should be enhanced relative to the other by the Lipkin mechanism [47] . 
It seems plausible that it should be C K +^ 1270 ^ because its main decay mode is via 
K + p° which can produce K +/ -y via vector meson dominance. 
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5 Pion polarizabilities and chiral symmetry con- 
straints 



5.1 Polarizabilities 

Polarizabilities are important observables which probe the structure of the pion. 
Using crossing symmetry and analyticity one can relate the amplitudes 77 — > 
7r + 7r~, 7r°7T at t — m\ and small s to the electric and magnetic polarizabilities 
of the charged and the neutral pion respectively. These pionic observables have 
been computed in Chiral perturbation theory at next to leading order [23| 121] and 
then at next-to-next to leading order [1H1 HHl [50j EI]- Some of these results which 
are relevant to our analysis will be recalled in the next subsection. There are simple 
relations between the subtraction constants introduced in the dispersive relations 
for the partial- waves and the electric/magnetic polarizabilities of the pion. From 
the partial waves with J = and J = 2 which are involved in our analysis we can 
access the dipole and the quadrupole polarizabilities. We list below the relevant 
formulas. The polarizabilities of the neutral pion are defined as follows from the 
expansions of H 1 ? , and if? around s — 0, t = m^ (e.g. [50J) 

2a H? , (s,t = ml) . n . s. n . 

s 12 
-2a H r '(s,t = ml) . n . s. n , 
m w s 12 

(the minus sign in from of H™_ is associated with our choice for the photon polar- 
ization vectors). The charged pion polarizabilities are defined by the expansion of 
the amplitudes after removing the Born term, i.e. defining H xx , = H xx , — H X y° rn , 

2a HlAs,t = ml) s. n . 

m n s 12 

-2a Hl(s,t = ml) n , s, n . 
m n s 12 

Performing the partial-wave expansion, the J = and J = 2 partial-waves are the 
only ones which contribute to the dipole and quadrupole polarizabilities. 

The polarizabilities are simply related to the subtraction constants in the disper- 
sive representations. Comparing with those fl23l) (1261) for Hq ++ and ++ we obtain 
for the polarizability differences 

(ai-^iW = -4 — (V2b^ + bV) 

(a 1 -p 1 ) lfl = ~ — (b^-^). (66) 
V3 m n 



20 



Similarly, the polarizability sums are obtained from the subtraction constants d 1 
which appear in the J = 2 spin-flip amplitudes h^, (see f[2"9"j) . f l3"Uj) ) 



(ai + /?i) 7r + = —10 am^ 

(a! + ft)^ = -10V2 am^d^ - v^rf^) . (67) 

Considering the quadrupole polarizabilities now, we list below the polarizabilities 
defined for 1 = and 1 = 2 amplitudes which can be easily combined using eq. (jl]). 
For the polarizability differences, we get 

(a a - fa® = — (Cl u 6C0) + firf + &'<*» + 10 m= C (°)) 

(a 2 - (3 2 )^ = ^ fog&C*) + 6'(2) + 10mM 2 )) . (68) 

For the sum of the quadrupole polarizabilities, we get 

(« 2 + /3 2 )(°) = 120vWa (6° d (0) + d'W) 

(a 2 + /3 2 ) (2) = 120^6^0 rf (2) + rf' (2) ) . (69) 

In these formulas the derivatives of the Omnes functions at s = appear, denoted 
by, e.g. (in. in e q- (I6"9"j) the quantities d'(°\ d'^ are given by sum rules 

W = L$_ + (70) 

where the integrals are defined in eq. ( 13"3"j) . Finally, the subtraction parameter 0$ 
which appears in the analysis through coupled channel unitarity is related to the 
1 = kaon polarizability as follows 

(Q1 - M f = ^L h f . (71) 
m K 



5.2 Constraints from chiral symmetry 

Chiral symmetry constrains the amplitudes 77 — > 7r°7r°, 7t + tt~ for small values of the 
Mandelstam variables s, t. Computations up to NNLO in the chiral expansion have 
been performed [HJ HH EH [51]. At this order, the amplitudes involve 13 coupling 
constants from the 0(p 6 ) chiral Lagrangian. More precisely, three combinations of 
such couplings are involved which we will denote d[, a 2 , b r for 7r°7r° and a[, a 2 , b r 
for 7r + 7r _ . Most of the 0(p 6 ) couplings are as yet undetermined, but in the case of 
7r°7T we can make use of a chiral sum rule for one coupling. We explain this below. 
We also recall the chiral expressions for the pion polarizabilities which allows one to 
assess which chiral constraints can be used in the fits to the experimental data. 
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a) 7r°7r°: 

The dipole polarizabilities of the neutral pion have the following expressions in 
ChPT, in terms of the coupling combinations a'[, a^, b r 



[Oil 



Pi) n 
Oil + 0i)tt o 



a 



167r 2 i* l2 m 7l 



1 m v 

"3 + 1^2 F 2 



86 r + X 1 _(^)) 



(72) 



(16tt 2 F2)2 

where the quantities Xi_(fi), involve chiral logarithms and known 0(p A ) 

couplings. Their detailed expressions can be found in ref . [50J . The chiral expressions 
for the quadrupole polarizabilities read 



(a 2 - (3 2 )n Q 

(a 2 + /^ttq 



a. 



a 



156 
45 



2b r )+X 2 ^)) 



5009 13435tt 2 16- 

I 1 h 

27 720 45 



(167r 2 F 2 ) 2 m w [ 27 ' 720 ' 45 " Z J ' ^ 

Based on resonance model estimates for the 0(p 6 ) couplings, numerical values for 
the polarizability differences at 0(p e ) were obtained [50]. We display them below 
together with the 0(p 4 ) values, which illustrates that 0(p 6 ) effects can be rather 
large for these observables: 



(«i - 0i)-K O 

(a 2 - (3 2 ) no 



0(p 4 ) 0(p 6 ) 

-1.0, -1.9 10 
20.7, 37.6 10 



" 4 fm 3 
" 4 fm 5 



(74) 



The reliability of naive resonance saturation models has not been established 
for 0(p 6 ) couplings. Here, we will only make use of a model independent estimate 
for the single coupling c r M . This estimate is based on a chiral sum rule [52j [53] 
associated with differences of correlators of two vector currents: (V 3 V 3 — V 8 V 8 ) or 
(yudydu _ yusysu^ f rom w hi c li the SU(3) coupling can be determined. Using 
three-flavour ChPT and matching to two-flavour ChPT it can be turned into an 
evaluation of the coupling Cg 4 (a simplified version of this sum rule was used earlier 
in ref. |54j). Such matching relations have been obtained recently by Gasser et 
al. [55] for the Lagrangian operators which do not vanish in the limit m u = ma = 0. 
We can obtain the matching relation for Cg 4 by going to this limit after taking a 
derivative with respect to m u + m^. Using the 0(p 6 ) ChPT calculations of the 
(V 3 V 3 ) correlator performed in ref. [56], and taking the derivative with respect to 
m 2 , one finds the following matching formula, 



"34 



192 x 167r 2 m 2 



61 



1 



K 



4 x 32tt 2 



iog^i + ; 



m + L r w ) 



+0(m 2 K ) 



(75) 
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(where m 2 K = lim r , 



=o itT'k)- This relation allows one to evaluate Cg 4 provided we 



further invoke a large N c argument which implies that Cg 2 (/i) should be suppressed 
compared with Cg-^/z) when /i ~ 1 GeV. The authors of ref. [53] have determined 
that 

n y3 (0) - II v k (0) = (1.92 ± 0.27) 10~ 2 (76) 

(using the notations of ref. [SB] for the (VV) correlation functions). Using further 
the explicit chiral formulas from this work [56J and the matching relation f[7"5l) above, 
we find 

c r M (fx = m p ) = (1.19 ±0.43) KT 5 . (77) 

Let us note the relation between Cg 4 and the combinations which appear in the 7r°7r° 
amplitudes 



8b r ) + 2(a r 9 - 2b r ) 



4094vr 4 c 



4 r 



34 



-34 



(A*)- 



(7f 



In other terms, the combination 6(ai — (3i) n o + m 2 (a 2 — Z^)^ depends on the single 
0{p & ) coupling, c 34 . Numerically, using \i = m p in the ChPT expressions, the 
following relation is obtained between the dipole and quadrupole polarizabilities 



6(«i 



with cJJ im 



f3 1 ) 7r o +m 2 r o(a 2 
: 4.75 ± 1.71. 



3: 



2)tt° 



(6.20 + 0.25. 



fiff, 
-34 1 



m. 



10" 4 fm 3 



(79) 
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b) 

In the case of the charged pion, the dipole polarizabilities are given in terms of 1$ — Z5 
at 0(p A ) and involve three combinations of couplings a[, a 2 , b r at 0(p 6 ) 



(ax - p x ) % + 

("1 + /3i)tt+ 



a 



3 l6 5J+ 16tt 2 F2 



86 r +X 1+ (^) 



0) 



(16tt 2 F2) 2 

The fully explicit expressions can be found in refs. [4*9| l5Tj. The quadrupole polar- 
izabilities are given as follows, 

.2 



(0L 2 ~ /3 2 )n+ 



a 



(a 2 + /3 2 )ttH 



lQT[ 2 F2:irfi 



a 



2 + 



l 6vr 2 F 2 



(12(2 



26 r )+X 2 _(/i) 



2062 10817tt 2 8 T 8 T 

1 1 h H 1 2 

27 1440 45 15 



Unfortunately, for the charged pion, there is no known model independent informa- 
tiorE on either of the three combinations a[, a r 2 -> b r . In this case, we will accept the 
estimates of ref. [51] stating that the polarizability difference should lie in the range 
(<*i - 0x)v+ G [4.70,6.70] 10" 4 fm 3 . 



4 A soft pion theorem due to Terazawa |57j has sometimes been applied to the 77 — > ir + ir~ 
amplitude. The theorem, however, applies to the amplitude 7* (9)7* (—q) —> 7r + (0)7r~(0) which is 
unrelated to 77 — > n + 7r~: it involves different chiral coupling constants. 
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6 Some details on the calculations 



6.1 Inputs for the tttt T-matrix 

We describe here our inputs for tttt scattering amplitudes which are needed for the 
S and the D waves. At medium energies tttt — > tttt phase-shifts and inelasticities 
have been measured in production experiments (see [58] for a review). Considerable 
progress has been achieved recently in measuring phase-shifts at very low energies 
as well as the I = 0,2 S'-wave scattering lengths PUJ dU [12] • In our analysis, 
the low-energy region is emphasized in the integrals by the s = subtractions. 
Interpolating between the medium and low energy regions is controlled by the set 
of Roy equations [59] . 
1=0 

Let us consider the S'-wave at first. As our main input, we will use the Roy equations 
results from ref. [HD] for the 1 = S-wsves in the energy range E < 0.8 GeV. In the 
range E > 0.8 GeV we perform fits to the experimental results of ref. |38^l . Below 1 
GeV, we will also make use of the phase-shift determinations of refs [63l 164"] . These 
also use the Roy equations as well as other dispersion relations as constraints. The 
phase-shifts differs from [BU] by a few degrees in the matching point region E ~ 0.8 
GeV: see fig [3J The differences in the results will serve us in estimating the errors. 

A well known feature of tttt scattering with / = J = is the sharp onset of 
inelasticity when the energy passes the KK threshold, which is caused by the /q(980) 
resonance. In the energy range which interest us here, it is a good approximation to 
ignore other inelastic channels and implement exact two-channel unitarity. We can 
deduce the required T-matrix elements directly from experimental inputs on tttt — > 
tttt and tttt — > KK scattering in the physical region and make use of analyticity in 
the unphysical region. Also, we work in the isospin limit and assume that this limit 
can be taken smoothly near the KK threshold. Then, the phase of the T-matrix 
element T 12 (s) = should be equal to the phase of Tn(s) = /o(s) (i- e - to the 

elastic tttt phase-shift) when yfs = 2m#. This threshold phase is actually not very 
well known at present: for illustration, the i^-matrix fit of Hyams et al. [38] gives 
5{2m K ) ~ 175° while the GKPY [61] analysis gives 5(2m K ) ~ 227°. Concerning T 12 , 
the two experiments by Cohen et al. [65J and Etkin et al. [66] are in disagreement 
close to the threshold, for both the phase and the modulus. The results of Cohen et 
al. are in accord with a large value of the threshold phase 5(2mjc) — 220° while, on 
the contrary, Etkin et al. favour a value smaller than 150°. An alternative possibility 
associated with the results of Etkin et al., is that of a fast variation of the phase 
in between the K°K° and the K + K~ thresholds, as suggested in ref. [67]. In that 
case, the isospin limit would not be smooth. 

5 It is likely that the phase-shift determinations in the region E > 1.5 GeV must be updated [txTl 
162) . This region plays a minor role in our analysis 
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Figure 3: S-wave I = tttt phase-shifts below 1 GeV: comparison of two 
determinations]^ [51]/ used in our analysis. 



In order to probe the sensitivity of the 77 data to the value of the threshold phase, 
we have performed our fits allowing it to vary in the range 150° < 5(2i71k) < 220°. 
Figure H] (upper plot) shows our fits of the tttt — > KK phase compared to the two 
data sets. Below the inelastic threshold, in the region 0.8 GeV < E < 2m^ we adopt 
the following simple description of the tttt phase-shift, involving a Breit-Wigner term 
plus a linear background 

Smiy/s) = a + bs + arctan "!° r o , 0.8 GeV < < 2m K . (82) 

m,Q — s 

Assuming given values for the phase-shift at E = 0.8 GeV and E = 2m k fixes the 
parameters a, b in terms of mo, Tq. These two parameters are then fitted to the 
experimental data in this region. A few values corresponding to different inputs for 
the threshold phase 5 7r7r (2m^) are collected in table SJ The mass and width of the 
/o(980) resonance in this simple parametrization are in reasonable agreement with 
the PDG values. The width is seen to be rather sensitive to the input threshold 
phase. The corresponding curves are shown on fig. HI At energies above 2mn, we 
describe both the tttt — > tttt and tttt — > KK phase-shifts by piecewise polynomial 
functions fitted to experiment. 

The modulus of T 12 (s) it needed in the unitarity equations (1T7|) in the unphysical 
region 4m^ < s < Am 2 K . Using analyticity and elastic unitarity it can be determined 
by the MO method. Since the left-hand cut of Ti 2 (s) can be expressed in terms of 
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Figure 4: Fits to the tttt — > KK phase (upper figure) and the tttt phase-shifts 
corresponding to different input values of the KK threshold phase. 



ttK — > ttK phase-shifts this MO equation is actually one component of the set of 
coupled Roy-Steiner equations (see e.g. [68]). We will employ here a simplified but 
reasonably accurate representation, 



T 12 (s) = (Aq + s{A + Bs + Cs 2 )) exp - / 

v ' 7T J Aral 



s r- 5 12 (s>) 



s'(s' - s] 



(83) 



where the two parameters A and Aq are chosen so as to reproduce the values 
T 12 (0) = 0.097 and its derivative T 12 (0) = 1.126 GeV^ 1 obtained from the full 
Roy-Steiner equations [68j and the two remaining parameters B and C are fitted to 
the experimental data in the range [1 — 1.5] GeV. 

The modulus of T 12 displays a peak associated with the / (980) resonance. As 
one can expect from the formula (183]) . the size of this peak is strongly correlated 
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S nn (2m K ) 



180° 200° 210° 220 



m (GeV) 
r (GeV) 



0.987 0.984 0.983 0.981 
0.056 0.039 0.033 0.028 



Table 4: Mass and width (in GeV) of the /o(980) resonance arising from fitting the 
parametrization / f&2|) to the data in the range 0.8 GeV < E < 2m k with different 
threshold phase inputs. 



with the value of the threshold phase Si^rrix) . This is illustrated in fig. The 
incompatibility between the results of Cohen et al. and of Etkin et al. below E ~ 1.2 
GeV is also apparent on this figure. Under the assumption of two channel unitarity, 
T12 is related to the inelasticity parameter rj^ in nir scattering by 



with o~p(s) — Jl — Amp/ 's. The results of Hyams et al. [38J on the tttt inelasticity 
which have relatively large error bars are compatible with both refs. (651 EE]. Once 
the 2x2 T-matrix is defined, the corresponding 2x2 Omnes matrix can be computed 
numerically (see [511 ES] f° r details). 

For the D-w&ve, we have relied on the two-channel i^-matrix representation of 
Hyams et al. [38] for representing the 7T7T phase-shift, updating the mass and width of 
the /2(1270) resonance to the PDG values. In this case the KK inelastic channel is 
not the physically dominant one below 1.3 GeV but is used as an effective description 
of inelasticity. Furthermore, the inelasticity quoted by the PDG at the energy of the 
/2(1270) is significantly smaller than the one obtained by Hyams et al. In practice, 
we have used a one-channel Omnes representation using the T-matrix phase instead 
of the S'-matrix phase in the inelastic region. 



In this case we have ignored inelasticity. For the S*-wave, we use the Roy parame- 
trization of ref. [60J below 0.8 GeV and make a simple fit to the data of refs. [HS1 ED] 
at higher energy. For the D wave we make a simple fit to the data found in the 
same references in the whole energy region. These fits are shown in fig. |6] 

6.2 Experimental data 

The main body of experimental data which we used are described in the publica- 
tions [HI [7] (charged pions) and [U |9] (neutral pions) by the Belle collaboration. 
They have measured differential cross-sections in the range cos# < 0.6 for 7r + 7r~ 




(84) 




1=2: 
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Figure 5: Modulus of the rnr — > KK S-wave T -matrix element using the represen- 
tation j53P fitted to the data, for several values of the KK threshold phase-shift. 



and cos 6 < 0.8 for 7r°7r°. We have taken into account all their data in the energy 
range E < 1.28 GeV i.e. 520 data points for 7r°7r° and 1152 data points for 7r + 7r~. In 
addition we have taken into account earlier experimental measurements of cross sec- 
tions integrated over cos 9 (in the same ranges as indicated above) from the Crystall 
Ball collaboration [7T] (7r°7r°) as well as Markll and Cello [72J [73] We have 

assigned equal weights to all the data points. Obviously then, the \ 2 is completely 
dominated by Belle's results. 

6.3 Parameters to be fitted 

The dispersive representations for the partial-wave amplitudes as written in sec. |3j 
involve 10 subtraction parameters. Not all of them will be determined from the fit. 

1) For the 1=2 D-waves we have actually assumed an unsubtracted dispersion rela- 
tion (i.e. the corresponding parameters and are determined from sum rules, 
see eq. (132])). 

2) We ave used some chiral constraints. Firstly, we have fixed the parameter b$ to 
be equal to its ChPT expression at one loop 

b m = JJ^+3A + 0(/) _ (85) 
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Figure 6: nn — > 7T7t isospin 2 phase-shifts: S-wave (upper curve) and D-wave (lower 
curve). 



Taking L$ and L\o from table 2 of ref. [74] gives 



&(°) ~ 



-0.40 ±0.30 



(86) 



As a further constraint we use the known value of the 0(p 6 ) coupling Cg 4 , as discussed 
in sec. I5.2I This provides one relation among the seven remaining parameters and 
leaves six parameters to be fitted. Finally, we perform the fits also imposing that 
the dipole polarizability difference of the charged pion lies in the range allowed by 
ChPT [51]. 



7 Results of the fit 

Differential cross-sections for 77 — > -k°7t° and 77 — > 7t + tt~ are evaluated with the 
J = and the J = 2 amplitudes computed as explained in sees. El HI EJ For 
77 — > 7t + tt~ all the J > 4 amplitudes corresponding to the Born term are also 
included. The values of the x 2 results after minimization, corresponding to the 
various data sets, and adding the statistical and the systematic errors in quadrature, 
are shown in table |5j These numbers correspond to a choice of left-cut cutoff A = —5 
GeV 2 and KK threshold phase 5{2rriK) = 200°. The values of the \ 2 f° r the various 
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Ndata 


X 2 /N d ata 


ttV Belle 

7r°7r° Crystal Ball 


520 
21 


1.26 
1.13 


7r + 7r~ Belle 
7r + 7T _ Cello 
n+ii- Mark II 


1152 
23 
67 


1.39 
1.23 
2.85 



Table 5: x 2 results from the constrained six parameters fit 



data sets are similar (which indicates compatibility) with the exception of the Mark 
II data [72J which shows some deviation. This feature was observed also in some 
previous analysis [75]. In more detail, the ability of our constrained dispersive 
representations to reproduce the experimental data is illustrated on fig. [7] (for 7r°7r°) 
and fig. [8] (for 7r + 7r~). The 7r°7r° amplitude is somewhat simpler than the ir + ir~ one 
due to the absence of the direct Born contribution in that case. This can be seen from 
the shapes of the differential cross-sections. In the energy region under consideration 
here only the S"-wave and two D-waves effectively contribute to the 77 — > 7r°7r° 
amplitude. In the energy region of the /2(1280) a significant contribution from the 
S'-wave background is needed for reproducing the experimental cross-section. Our S- 
wave amplitude based on two-channel unitarity cannot be quantitatively trusted in 
this region. Figure [9] illustrates the role of the various partial- waves in the integrated 
cross section. 







°K 


6 '(2) 


c° 


d° 


Fit 


-7.99 


-7.16 


3.47 


0.22 


-0.48 


SR (A = -5) 
SR (A = -3) 


-6.00 
-7.62 


-5.61 
-6.84 


3.29 
3.22 


0.19 
0.10 


-0.54 
-0.57 



Table 6: Comparison of the subtraction parameters as determined from a fit of the 
experimental data and as determined from sum rules, as a function of the cutoff A 
on the left-cut integration. Units are in appropriate powers of GeV. 



The dispersive representations are based on over-subtracted dispersion relations. 
Assuming reasonable high-energy behaviour, five of the subtraction parameters 
could be written as sum rules (see sec. 13. ip . In practice, the result of such sum 
rules depend on the cutoff A that one introduces on the left-hand cut integration 
(since the precise behaviour of the integrand for large negative s is not known). 
Table [6] shows the numerical (central) values of these parameters, as generated by 
the fit, and the sum rule evaluations. The table shows that, for physically reason- 
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E (GeV) 



Figure 7: Differential cross-sections (systematic and statistical errors added in 
quadrature) and integrated cross-section (Belle's systematic errors are shown sepa- 
rately) for 77 — > 7r°7T . The solid line is the result of our constrained fit. 
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Figure 8: Same as fig. [7| for 77 — > 7T+7T 
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Figure 9: Contributions from the S and the D partial-waves to the 77 — > 7r°7r° 
integrated cross-section 



able values of the cutoff, A ~ —3, —5 GeV 2 , the fitted values are qualitatively in 
agreement with the sum rule ones. 

• /o(980) region: 

Our parametrization of the input nn phase-shifts allows for some freedom to vary the 
value of the phase-shift at the KK threshold. While the overall x 2 is hardly sensitive 
to this small energy region, Belle's data provides a detailed picture of the / (980) 
peak because of the large statistics. A comparison of our results corresponding to 
different values of the threshold phase-shift <5(2m^), with Belle's data, is illustrated 
on fig. [IOJ For clarity of the figure, the systematic errors are not shown, in this 
region they are of the order of 12 nb for 7r + 7r~ and 3.5 nb for 7r°7r°. In the case 
of 7T 7r°, comparison with the data favours values of the threshold phase <5(2mx) ^ 
180°: for smaller values the peak is too flat and displaced to the right. The value 
<5(2mx) — 200° eventually provides the closest agreement with the experimental 
peak. The shape of the /o(980) in 7r°7r° and in 7r + 7r~ is predicted to differ because 
of the different sign of the interference between the 1 = resonant amplitude and 
1 = 2 amplitude. There is some indication of this feature in the data. Comparing 
the MO results with the central values of the data for 7r + 7r~ one must keep in mind 
that the systematics are larger than for 7r°7r°. However, essentially the same value 
of the threshold phase S(2rriK) — 200° also provides the best agreement with the 
shape of the / (980) peak. When 5{2m K ) gets smaller than ~ 170° the structure 
resembles a cusp rather than a peak. Belle's statistics are very high for 7r + 7r _ and the 
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Figure 10: Integrated cross-sections in the region of the / (980) peak. The solid, 
dashed and dotted curves are the result of our MO representation corresponding to 
three values of the tttt phase-shift at the KK threshold. The long-dashed curves 
correspond to using the phase-shifts of ref. The experimental results of Belle 
\EI are shown with their statistical errors only. 



energy bins AE = 5 MeV are smaller than the spacing between the K°K° and the 
K + K~ thresholds ( ~ 8 MeV, the two thresholds are indicated by arrows in fig. fTOi) . 
Belle's data are compatible with a small isospin breaking at the KK threshold since 
agreement with the MO amplitude (which has no isospin breaking) in this region 
is not worse than elsewhere. A strong isospin breaking scenario has been proposed 
by Au et al. [57], according to which the tttt — > KK phase drops very sharply in 
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between the two thresholds. Such a scenario is also not ruled out by Belle's data. 
Isospin breaking near the KK threshold involves ao(980) — /o(980) mixing. Studies 
of how this affects the shape of the ao(980) peak have been performed (e.g. [76], [77]). 
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Figure 11: Integrated 77 cross- sections at low energy. The experimental results 
from [711 172} are shown and also the results from ChPT calculations as given in 
/5U[ 151} . The solid line is the result from the MO amplitudes. 



• Low energy region, pion polarizabilities 

Next, we consider the low energy region. Fig. [TT1 shows our result for the integrated 
cross-sections in this region and also shows, for comparison, the result from ChPT 
calculations at NLO and NNLO. As another comparison, we show in fig. [T2] the 
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Figure 12: Amplitude H ++ (s,t = u)/s in the sub-threshold region and \H ++ (s,t = 
u)\/s above the threshold (as in fig. 8 of ref. J48J ) compared to ChPT results. 



amplitude H ++ (s)/s in the sub-threshold region normalized as in fig. 8 of ref. jlH] 
(above the threshold the modulus is shown). In this region, the dispersive amplitude 
lies rather close to the chiral amplitudes. In particular, it has an Adler zero close to 
s = as has been anticipated in ref. [21] ■ The differences between the dispersive 
results and ChPT are to be attributed to effects of chiral order 0(p s ) and also to 
some differences in 0(p 6 ) coupling constants (see below). 

Expanding the 77 amplitudes around s = with t = m^, one accesses pion 
polarizabilities (see sec. [5]). The results for these quantities deduced from our fitted 
amplitudes are collected in table [7J and compared with the results from ChPT at 
NNLO as given in [5DJ EI]- Our fit was performed with the constraint that for 
the charged pion, the polarizability difference (cti — f3i) n + should lie within the 
range of the ChPT calculation. The fit prefers the lowest value in the allowed 
range. This trend differs from the result of Fil'kov and Kashevarov [79] who fitted 
the charged amplitude only. For the neutral pion, we have imposed a constraint 
between the dipole and the quadrupole polarizabilities (see (1791 ). The result for 
the polarizability difference (cui — /3i) n o is then in acceptable agreement with the 
ChPT prediction. Our results for the dipole polarizability sums are in agreement 
with ChPT. For the quadrupole polarizability differences, our results are slightly 
smaller than ChPT. This corresponds to somewhat different results for the chiral 
coupling constant combinations d[, a r 2 and h r for which simple models have been used 
in ref. [501 loTj . Table [S] shows the values of these constants resulting from the fit and 
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7T U 


Fit 


ChPT ±Res. mod. 


(«! - (3i) n o 

(a 2 - f3 2 ) n o 
(a 2 + (3 2 )nO 


-1.25 ±0.08 ±0.15 
i oo _i_ n i o _i_ n no 

32.1 ± 0.9 ± 1.9 
-0.19 ±0.02 ±0.01 


-1.9 ±0.2 

1 1 ± n 3 

37.6 ±3.3 
0.04 


7T+ 






(«i - /3i)tt+ 

(«1 + /?l) 7r + 

(«2 - fti)^ 
("2 + 


4.7 

0.19 ±0.09 ±0.03 
14.7 ± 1.5 ± 1.4 
0.11 ±0.03 ±0.01 


5.7 ± 1.0 
0.16[0.16] 
16.2(21.6] 
-0.001 



Table 7: Results for dipole and quadrupole polarizabilities (in units of 10 -4 fm 3 and 
10 -4 fm 5 respectively) of the ir° and ir + compared with the values from ChPT at 
0(p 6 ) associated with a model for the LEC's (from ref. [50]/ for the ir°, and fBl] for 
the 7r + . The numbers in brackets correspond to using the ENJL model \ 78} for the 
LEC's). The central values shown correspond to a fit with a left integration cutoff 
A = —5 GeV 2 . The first error corresponds to the variation of the input parameters 
and the second error reflects the uncertainties in the experimental 77 data. 



compared with those from a resonance model and also from the ENJL model [78J. 
We note that, in the NJL model, the dipole and quadrupole polarizabilities can be 
calculated directly, e.g. [HTil EI]- Our predictions differ from 0(p 6 ) ChPT for the 
sums of the quadrupole polarizabilities. The 0(p 6 ) chiral couplings cancel out in 
the expressions of these observables. As a consequence, they are expected to be 
sensitive to effects of chiral order p 8 [51] . 

Let us also mention that our results for (a 2 — (3 2 ) n o tn + are somewhat smaller than 
those obtained in refs. [821 [79] , 

(« 2 - AO.o = (39.70 ± 0.02) 10" 4 (a 2 - f3 2 )^+ = (25.0+J^) 10~ 4 (87) 

(in units of fm 5 ) from fitting subtracted dispersive representations (they used a 
combination of s- fixed and t- fixed dispersion relations) of 77 — > 7t + tt~ (ref. |79j). 
77 —?• 7r°7T (ref. [82]), amplitudes to experimental data. This difference in the results 
is to be attributed, we believe, to our combining 7r°7r° and 7i + ir~ data in the fits as 
well as our using a more sophisticated treatment of the final-state interaction in the 
S'-wave, which plays a crucial role for polarizability differences. Our results for the 
quadrupole polarizability sums (a 2 + (3 2 ) n o^+, which are controlled by the D-waves, 
are in rather good agreement with [821 [79] . 
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7T° 


a\ a T 2 b r 


Fit 

Res. mod. [ENJL] 


-25.9 ± 1.6 ±3.7 8.6 ±0.8 ±1.8 3.4 ±0.4 ±0.1 
-39 ±4 [-23.3] 13 ±2 [14.9] 3 ±0.5 [1.7] 


7T+ 


a r 2 b r 


Fit 

Res. mod. [ENJL] 


-25.0 ±2.2 1.4 ±1.8 ±1.4 0.2 ±0.3 ±0.1 
-3.2[-8.7] 0.7[5.9] 0.4[0.4] 



Table 8: Values of the three combinations of 0(p 6 ) coupling constants at the scale 
p = 0.77 GeV involved in the 77 — > 7r°7r° amplitude (second and third row) and 
in the 77 — > tt + it~ amplitude (fourth and fifth row). The values from a resonance 
model used in refs. f5U\ 151} and from the ENJL model J 78} are compared with the 
values deduced from our fitted dispersive amplitudes. Errors are as in table [71 



• Errors 

The errors quoted in tables [71151 have been estimated as follows. The uncertainties 
associated with the description of the left-hand cut has been evaluated by varying 
all the input coupling constants as well as the integration cutoff A which was varied 
between —3 and —10 GeV 2 . Concerning the final-state interaction, we have varied 
the vrvr scattering lengths (ag = 0.220±0.005, a 2 , = -0.0444±0.0010) and used, below 
E = 1 GeV two different representations for the 1 = S-wave tctt phase-shift [60, 64J. 
We have also varied the resonance parameters in the D-wave. Concerning the errors 
in the fitted parameters, the usual criterion based on increasing the \ 2 by one unit 
is based on the assumption that the experimental errors are statistical and that the 
correlation matrix is known. These assumptions are not valid in the present case. 
We have therefore adopted a more phenomenological criterion, considering the x 2 
per point instead of the total one and allowing it to increase by 0.5. The errors 
associated with the fitted parameters and those associated with the input data are 
quoted separately in tables ITfHl The ir + polarizability difference (ai — f3i) n + is not 
completely determined by the fit as it lies at the boundary of the allowed value. No 
error can be quoted in this case. Correspondingly, we quote a single error for the 
LEC combination a\ in table [SJ 

8 Conclusions 

We have reconsidered the MO dispersive representations of photon-photon scat- 
tering amplitudes and applied them to the recent results of the Belle collabora- 
tion |1 [S Ij. This method is general and follows from the non-perturbative 
features of QCD. Its range of applicability can be extended up to slightly above 1 
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GeV by taking into account the main source of inelasticity in irir scattering. Our 
description of the left-hand cut includes the contributions of the vector, axial- vector 
as well as tensor resonances i.e. all relevant resonances with masses up to 1.3 GeV. 
We found that the tensor resonances play a significant role in the left-hand cut. We 
employ MO representations somewhat different from previous works, keeping the 
left-cut parts in the form of spectral integrals. In this manner, one avoids polyno- 
mial ambiguities associated with propagators of particles with non-zero spin. Fur- 
thermore, the spectral function displays in a clear way how a regularization occurs 
in exchanges of resonances of different types. 

We have also argued that chiral constraints can be imposed on the subtraction 
constants using model independent information on p 6 chiral couplings. We have 
shown that a sum rule on a single p 6 parameter for 7r°7r° provides a non trivial con- 
straint. Making use of all these theoretical constraints, one obtains an amplitude 
containing only six parameters. We introduce a cutoff on the left-cut spectral inte- 
gration, accounting for a Regge-type regularization, but no other cutoff. In contrast 
to ref. [27], the QED Born term is used unmodified by form-factors. 

We have fitted the subtraction parameters to a data set containing 541 points 
on 77 — > 7r°7r° and 1242 points on 77 — > 7r + 7r~. This set is largely dominated by 
Belle's results. We find good compatibility between Belle data on 7r°7r° and the 
lower energy data from the Crystal Ball collaboration. In the case of 7r + 7r~ we find 
compatibility between the data from Belle and the data from Cello as well as the 
data from Markll below 0.4 GeV. The compatibility with the Markll data in the 
range 0.4 — 1.0 GeV is more marginal. 

In the region of the KK threshold and the /o(980) peak, we find reasonable 
agreement between Belle's data and our parametrization which assumes a smooth 
isospin limit. The charged channel, however, does not rule out the possibility of 
some isospin violating effects in the shape of the peak. The cross-sections around 
the peak in both the neutral and charged channels are best reproduced for values 
of the threshold phase 5(2mx) — 200° ± 20°. From our calculation we also obtain 
the amplitude 77 — > KK for / = J = 0. It is unfortunately difficult to probe this 
part against the experimental data on 77 — > K + K~ , K°K° since the experimental 
amplitudes contain admixtures from isospin 1 = 1. 

We have shown that Belle's data are compatible with the pion polarizabilities 
predicted in ChPT. In the case of the 7r° we have derived a refined value. We 
have also derived results for the quadrupole polarizabilities of both the charged 
and neutral pion. Our results for the differences «2 — 02 are somewhat smaller than 
those derived in refs. [821 [79] . We believe this to be due to our using more precise nn 
phase-shifts at low energies for the S"-waves. It is clear, however, that experimental 
data on photon-photon scattering at low energies (i.e. below 0.5 GeV) are most 
efficient for determining pion polarizabilities with precision. It is hoped that such 
experiments will be performed at facilities like KLOE2 or BESIII. 
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